Highly anisotropic superconducting gap in underdoped Bai_ a; K :! ;Fe2As2 
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The in-plane London penetration depth, AA(T), was measured using a tunnel diode resonator 
in single crystals of Bai- :E K :E Fe2As2 for five doping levels x ranging from underdoped, T c =ll K 
(a;=0.17), to optimally doped, T c =38 K (a? pt= 0.35). In the optimally doped samples, AA(T) shows 
exponential saturation in T —¥ limit suggesting a fully gapped superconductivity. The lowest-Tc 
samples show much stronger non-exponential variation of AA(T). Fitting the data to a power-law, 
AA(T) = AT n , reveals a monotonic decrease of the exponent n with x towards the underdoped 
edge of the superconducting dome. Comparison with n w 1.2 reported in KFe2As2 (T c =3.5 K , 
x = 1), suggests a dome-like variation of n with x, implying an evolution of the topology of the 
superconducting gap from full and isotropic in the center of the dome towards strongly anisotropic 
and, eventually, nodal at the dome edges. 



The experimental determination of the symmetry of 
the superconducting gap is important for understanding 
the mechanism of superconductivity in iron-based super- 
conductors 0, Q • Measurements of London penetration 
depth , thermal conductivity @, 0] and specific heat 
[ItKj] in electron doped Ba(Fei_ a; Co a: )2As2 (BaCol22) 
suggest that superconducting gap shows strong evolu- 
tion with doping, developing nodes at the dome edges 
S E > 12 1 • This doping-evolution is consistent with ob- 
servations of a fully gapped superconductivity in effec- 
tively close to optimally-doped LiFeAs fl3l - [l6j ] and nodal 
superconductivity in effectively overdoped KFe2As2 (l7r - 
Il9j . It is also consistent with predicted doping-evolution 
for the s± model 12, 21|. On the other hand, nodal be- 



havior is observed at all doping levels in isovalently sub- 
stituted BaFe 2 (Asi_ a; P a: )2 (BaP122) [5(J. This notewor- 
thy difference in two systems based on the same parent 
compound prompts a detailed study of the hole doped 
Bai_ a: K :c Fe2As2 (BaK122). The superconducting gap in 
BaK122_was studied intensivel^usdng ARPES |22l. l23jj. 



NMR [24| . penetration depth [25j, |26j and thermal con- 
ductivity 0, l27j |. however, no systematic doping - de- 
pendent study reaching the dome edges was undertaken 
so far. 

In this work we study the evolution of the tempera- 
ture dependence of in-plane London penetration depth, 
AA(T), in high quality single crystals of Bai_ a: K 2 ;Fe2As2. 
We find that the optimally doped samples show expo- 
nentially weak temperature dependence in T — > limit, 
suggesting a fully gapped superconductivity. This con- 
clusion is consistent with the temperature-dependent su- 
perfluid density in these samples, which can be well fitted 
using self-consistent 7-model with two full gaps in the 
clean limit [28]. The lowest-T c samples show an excep- 
tionally strong sub-quadratic temperature dependence. 
Fitting the experimental AA(T) below T c /3 to a power- 
law, AA(T) = AT n , we find a monotonic decrease of the 



exponent n with concomitant sharp increase of the pre- 
factor A towards the low x edge of the superconducting 
dome. Comparison with close to T-linear behavior found 
in heavily overdoped KFe2As2 [If|, suggests a universal 
development of nodes at the edges of the superconduct- 
ing dome in both electron- and hole-doped BaFe2As2 - 
based superconductors. 

Single crystals of Bai_ :c K :E Fe2As2 were grown using 
high temperature FeAs flux method [29|. AA(T) was 
measured using tunnel-diode resonator technique |3oLl3lj . 
Placing a sample into the inductor causes the shift of 
the resonant frequency, A/(T) = — G4nx(T). Here 
4-7rx(T) is magnetic susceptibility and G is a calibra- 
tion constant determined by physically pulling the sam- 
ple out of the coil. With the characteristic sample size, 
R, Attx = (A/_R) tanh(i?/A) — 1, from which AA can be 
obtained [3l|, [32|. The excitation field in the inductor, 
H ac 20 mOe, is much smaller than H c \. 

To compare sharpness of the superconducting tran- 
sition, Fig. [Ha) shows normalized RF susceptibility of 
Bai_ a; K a ;Fe2As2 samples used in this study. The super- 
conducting transition remains quite sharp even for the 
most underdoped samples where T c (x) is very sensitive 
to small variations of x. The values of x were determined 
by the empirical fit [33| of the experimental T c {x) data 
[29 . l34| . The values of T c were determined from the posi- 
tion of the maximum in the first derivative, dAX(T)/dT. 
In our samples we obtained T c of 11.2, 14.5, 18.6, 30.0, 
and 38.7 K, corresponding to potassium concentrations 
of a; =0.17, 0.18, 0.20 (all ±0.01), 0.28 (±0.02), and 0.35 
(±0.03), respectively. The low-temperature variation of 
AA(T) up to T c /3 is shown in Figs. \Eb) and (c). Fig- 
ure [UJb) compares the data for limiting compositions 
x=0.17 and 0.35, revealing a big difference in the magni- 
tude of AA(T) . Two curves for pure KFe2As2 are shown 
for reference [19(. Figure QJc) shows AA(T) on the same 
scale for all concentrations. The data are offset for clar- 
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ity and red lines represent the power-law fit. However, 
a closer look shows significant deviations of the data in 
heavily underdoped (x — 0.17) and in optimally doped 
(x = 0.35) samples. At lowest temperatures, it becomes 
significantly sub-quadratic for the former and closer to 
exponential for the latter. 

First we attempted to fit the data for two high- 
est T c samples to the single gap s-wave BCS function, 
AA(T)/A(0) = ^7rA /2fc s T exp(- A /fc B T), where A 
is the size of gap at T — 0. The Ao values from the best 
fittings are 0.73 k B T c and 0.87 k B T c for x = 0.28, and 
0.35, respectively. While the fit quality was good, both 
Ao values are much smaller than in single full-gap super- 
conductors where Ao = 1.76 ksT c . Such small gaps are 
expected in superconductors with A m j n < A max , cither 
due to gap angular variation (anisotropy) or variation 
between different Fermi surface sheets. 

A standard way to analyze AA(T) is to fit it from the 
lowest temperature up to T up » T c /3. In a single -gap 
s-wave superconductor this limit is determined by reach- 
ing nearly constant value of the superconducting gap A , 
below which the temperature dependence is exponential. 
For various nodal gaps, the dependence is expected to be 
power-law, T-linear for line nodes and T 2 for point nodes 
in clean limits. For the anisotropic gap or multi-gap su- 
perconductors with the variation of the gap magnitude 
over the Fermi surface between A max and A m j n , the T up 
is determined by A m j n , while T c by A max , so that T up 
range of characteristic temperature dependence can be 
smaller than T c /3. 

We therefore check the alteration of the fitting param- 
eters by choosing different temperatures for the upper 
limit, T up < T c /3. The dependence of n and A on T up is 
shown in Fig. [2][b). The highest-and lowest- T c samples 
exhibit monotonic increase and decrease of n on T up — > 0, 
approaching very different values of 4 and 1.5, respec- 
tively. 

The exponents n for samples with x = 0.18, 0.20, and 
0.28 do not show a significant variation with T up < T c /3, 
indicating robust power-law behavior, but show system- 
atic increase of n with x. The decrease of n with decrease 
of x can be clearly seen by in the top panel of Fig. [2j in 
which all data are plotted vs (T/T c ) 2 - 78 , where n = 2.78 
is the exponent for the optimally doped samples. The 
dependence of the power-law pre-factor A on T c was an- 
alyzed by fixing n=2.0, 2.3, and 2.5 and is shown in the 
inset revealing a significant increase with decreasing T c . 
The A value for sample x=0.17 is 30 nm/K 178 , out of 
scale for the plot and is not shown. 

A smaller than weak - coupling value of A m i n obtained 
from low-temperature BCS formula implies two - gap su- 
perconductivity and the analysis must be extended to 
the full - temperature range. The most convenient quan- 
tity is the supcrfluid density, which can be calculated 
from the first principles. In the optimally doped samples, 
we fit the data using clean-limit 7— model [28j. Symbols 
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FIG. 1. (a) Normalized AA(T). (b) and (c) Low-temperature 
parts of AA(T). Experimental data are displayed with solid 
circles. Solid red line represent power-law fitting curves. Two 
dashed lines are data for pure KFe2 As2 from Ref. [l9|. 



in Fig. H show superfluid density, p„(T) = A 2 (0)/A 2 (T), 
for the sample with x = 0.35 calculated from A(T) = 
AA(T) + A(0) with A(0) = 200 nm [Hj]. Solid lines show 
self-consistent 7-model fit for two-full-gap superconduct- 
ing state [2^| with p s — jpi + (1 — 7)^2, where p\ and pi 
are partial superfluid densities. Insert shows two super- 
conducting gaps Ai and A2 calculated during the fitting 
procedure. The estimated gap values are 6.5 and 3.3 
meV. Specific heat jump produced the value of ~ 6 meV 
for the larger gap 36] . 

Upon departure from optimal doping, the exponent 
n shows notable evolution with x decreasing from 4 to 
about 1.5. London penetration depth is very sensitive 
to pair-breaking disorder, modifying AA(T) at low tem- 



peratures [28|, |37J. In Bal22 - derived compounds it was 
also suggested experimentally (2|| (2|| [H, |3{|. Within 
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FIG. 2. (a) AA vs. (T/T c ) 2 - 8 up to T c /3. (b) Variation of ex- 
ponent n obtained from a power-law fitting, AA = A(T/T c ) n , 
as a function of the upper end temperature of the fitting 
range. Inset: A vs. T c . The three curves represented by 
triangle, circle, and square were obtained with hxed n =2.0, 
2.3, and 2.5, respectively. 




FIG. 3. Symbols: superfluid density p a (T) = A 2 (0)/A 2 (T) 
calculated with A(0) = 200 nm in Bao.65Ko.3sFe2As2 [3a |. 
Solid lines represent the ht to a two-gap 7 model, p s = ypi + 
(1 — r y)p2- Inset: superconducting gaps Ai(T) and Ai(T) 
calculated self-consistently during the fitting. 



clean limit 




A(T) should be exponential in the 
However, pair-breaking scattering 
(which in this case can be caused by non-magnetic impu- 
rities and dopant ions) turn the behavior into a power- 
law with the exponent n approaching 2 in the dirty limit 
37 . l42l [43| ] . Since the superconductivity in BaK122 is 



induced by doping, we cannot ignore the effect of dis- 
order on the variation of exponents. However, it would 
be natural to expect increase of scattering with x, and 
thus decrease of the exponent, opposite to the trend in 
our data. Similarly, disorder effect cannot explain nodal 
state in the end member of BaK122, very pure KFe2As2 
with n=1.2 [3, El- In addition, our most underdoped 
sample shows the exponent n = 1.5 clearly well below 
the limiting value of 2 for pair-breaking scattering. Thus 
we conclude that the variation of the exponent n, found 
in our study, is caused by the changes in the supercon- 
ducting gap structure with doping. 

The evolution of the power-law behavior in 
Bai_ x K 2; Fe2As2 superconductors is summarized in 
Fig. IHa). Solid circles show exponent n with the 
error bars estimated from the fitting to the different 
temperature ranges (such as shown in Fig. [U(b)) and 
open circles show the pre-factor A calculated for a fixed 
exponent n — 2.3. Also shown are the exponents for two 
stoichiometric (clean) compounds, KFe2As2 [l3| and 
LiFeAs [15j. The dashed line represents our picture of 
the exponent variation with doping that, in our opinion, 
reflects developing anisotropy of the superconducting 
gap. To relate to the phase diagram, Fig. |4jb) shows 
magnetic and superconducting transitions vs doping 
from neutron scattering [HI]. The pre- factors of the 
power-law fit show a sharp increase in the antiferro- 
magnetic (AFM) region, similar to FeCol22 [3], which 
indicates microscopic coexistence of superconductivity 
and long-range magnetic order [4o| . 

The value of the exponent n—1.5 in the sample with 
:r=0.17 suggests the existence of line nodes in the su- 
perconducting gap in heavily underdoped compositions. 
This value is too small to be explained in the full-gap 
s± scenario with strong pair breaking scattering (38|. In 
superconductors with line nodes the exponent n varies 
between n=l in the clean limit and n=2 in the dirty 
limit [46| . Our values of n are well within this range and 
outside the dirty s±. 

We would like to emphasize that in electron-doped 
Bal22 a clear signature of nodes appears only for the 
c-axis as shown by the penetration depth [U and thermal 
conductivity measurements 0]- Here we find that in the 
hole-doped Bal22 the nodes seem to appear even in the 
a&-plane measurements. This may be indicative of the 
different nodal structure in these two systems. In both 
cases, the nodes develop far from the optimal doping. 

Strong dependence of the superconducting gap struc- 
ture on the shape of the three-dimensional Fermi surface 
is predicted in the s± model 11(. In the same model, 
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FIG. 4. (a) n vs. x diagram. The red solid circles represent n 
from various temperature range between T c /6 and T c /3. The 
exponent n ~ 1 and 3.1 in KFe2As2 [hJ and LiFeAs [l5l |. 
respectively, are from elsewhere, (b) T-x phase diagram. The 
AFM region is from neutron scattering data [33]. 



the gap also evolves with doping due to changing nest- 
ing conditions, in rough agreement with the evolution 
we see here 



21| . However, the nodal superconducting 



gap can also be explained by the phase transition from 
fully-gapped s± to nodal <i-wave, which have relatively 
small difference in ground state energy [21|. However, 
this model does not explain a gradual evolution of the 
exponent n with doping. 

In conclusion, the measurements of A(T) in 
Bai_a;K x Fe2As2 suggest doping-dependent evolu- 
tion of the superconducting gap from full and isotropic 
at optimal doping to highly anisotropic and, eventually, 
nodal in heavily underdoped samples. 
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